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If a tennis ball is held above a basket ball with their centers vertically aligned, and the balls are 
released to collide with the floor, the tennis ball may rebound at a surprisingly high speed. We show 
in this article that the simple textbook explanation of this effect is an oversimplification, even for 
the limit of perfectly elastic particles. Instead, there may occur a rather complex scenario including 
multiple collisions which may lead to a very different final velocity as compared with the velocity 
resulting from the oversimplified model. 

PACS numbers: 45.70.-n,45.70.Qj,47.20.-k 



I. INTRODUCTION 

Consider a set of two balls made of the same viscoelas- 
tic material whose centers are vertically aligned at posi- 
tions z\ and Z2 as sketched in Fig. [T] Let R\ and i?2 
be the radii of the particles and A/i their initial vertical 
spacing. At time t = we release the particles to collide 
with the floor. The question for the maximum height 
reached by the upper sphere after the collision is then 
a common textbook problem. The textbook solution is 
based on the assumption that the collisions of the lower 
particle with the floor and the subsequent collision of the 
lower particle with the upper one are separate two-body 
interactions which may be treated independently, that 
is, one disregards the duration of the collisions. Whether 
this assumption is justified depends, of course, on the 
initial conditions, the particle sizes and the material pa- 
rameters. 

The experimental investigation of the described prob- 
lem is tricky as even a small deviation from the vertical 
alignment of the initial positions of the spheres leads to 
considerable post-collisional velocities in horizontal direc- 
tion, in particular, if more than two balls are involved. 
Simple but effective techniques were introduced [TH3], 
which allow for experiments with chains of up to about 




FIG. 1: Sketch of the problem. 



7 aligned balls. 

In this paper we will show that although the problem 
looks simple, there may emerge rather complex behavior 
including multiple collisions. The simple solution men- 
tioned above is, thus, only a special case of a more general 
solution. 



II. INDEPENDENT COLLISION MODEL 
A. Collision Scenario 

To introduce our notation and for later reference let us 
first reproduce the solution for the independent collisions 
model (ICM), that is, the final velocity of the upper ball is 
the result of a) the inelastic collision of the lower ball with 
the floor and b) the inelastic collision of the lower ball 
with the upper ball. Note that here and in the following 
we restrict ourselves to the situation where mi > m^- 
For the opposite case, m\ < mi, it may be shown that 
the sketched sequence of collisions fails, even under the 
assumption of perfectly elastic collisions [TJ 0]. 

Inelasticity of the balls is described by the coefficient 
of restitution which relates the precollisional relative ve- 
locity Vij = Vi — Vj of colliding particles, i and j to the 
post-collisional one, v'^ — v[ — u'-, 

e = -v'ij/ v i:j . (1) 

The lower sphere (mi) reaches the floor at time iio = 

^- g {zf ] at the velocity v w = -^2g{zf ) - R x ) 

from where it is reflected with v[ = —eviQ. The first 
index always denotes the considered particle and the sec- 
ond its collision-partner. Index stands for the floor, 
1 for the lower and 2 for the upper sphere. Upper in- 
dex (°) stands for initial values. Post-collisional values 
are marked by primes. The particles collide then at 
ti2 = tw — t, 10 fiVe) = + ^ whereby the upper sphere 
is located at position z 2 i = z 2i = — ^#^12 + z 2°' > with the 
initial height = z^ +Ri +R2 + Ah. At £12 the lower 
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particle has the velocity vi 2 — —gAt + v^Q and the upper 
f 2 i = —gt\2- Employing the collision rule, Eq. ([!]), and 
conservation of momentum, we obtain the final velocities 



m 2 V2i + roikvi2 + e(v 12 - v 2 i)] 

mi + m 2 
m\v 12 + m 2 [v 2 i + e(v 2 i - vi 2 )j 



mi + m 2 



and the relative velocity 
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In the case of elastic (e — 1) and instantaneous (Ah — > 0) 
collisions we find 
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For fj, = 1/3 the lower sphere loses all its kinetic energy 
(v[ = 0) and the upper sphere rebounds with twice its 
initial velocity. In the limit \x — > we recover the well 



known textbook result 



-3«io, that is, the upper 



ball rises to about nine times of its initial dropheight 
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9z 2 0) ~8(2R 1 +R 2 )). 



The system of bouncing balls can be exhaustively de- 
scribed also for more than two spheres [S], provided the 
collisions are considered as isolated events, that is, only 
two-particle interactions are taken into account. 



B. Coefficient of Restitution Resulting from the 
Solution of Newton's Equation 

The contact of viscoelastic spheres is described by the 
(modified) Hertz contact law [5]. In this article we will 
use a simplified force law since it allows for a exhaus- 
tive analytical solution of the problem. To justify this 
approximation, we will show later by means of numeri- 
cal simulations that the more correct Hertz contact force 
leads to qualitatively identical results, see Appendix |B| 

We describe the contact of dissipatively interacting 
particles, i and j, by 



F 



(7) 



as a function of the mutual compression 



elastic force which would lead to a resulting erroneous at- 
tractive force, see e.g. [7]. Therefore, the min[. . .] func- 
tion is applied to take into account that the interaction 
force is always repulsive (negative). 

Consider an isolated pair of colliding particles i and j 
approaching one another at impact rate v = £(i = 0) at 
t = 0. Using the force, Eq. Q, we obtain the relative 
velocity after a collision by solving Newton's equation of 
motion, 
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ij SiJ 



(9) 



with the effective mass mfj = miirij/ (mi + rnj) and ini- 
tial conditions £ii(0) = and £ij(0) = v. The collision is 
complete at time t c when £y(i c ) = [7]. 

Of course, for a pairwise collision the final velocity as 
obtained from Eq. ([I]) must coincide with the final ve- 
locity as obtained from integrating Newton's equation of 
motion. Therefore, the solution £ij(t c ) of Eq. ^ allows 
to relate the coefficient of restitution e to the parameters 
k and 7 of the force law, Eq. ([7|, via 



e = — 



(10) 



Straightforward calculation '7\ yields for the duration of 
the collision 
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For the coefficient of restitution we obtain 
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where D, = yJP 



(13) 

Note that e depends on the pa- 
rameters of the force and the effective mass m of the 
colliding particles, that is, e = s(k,j,m cB ). Thus, s may 
not be considered as a pure material constant. 



= max [0, R, + Rj - \f t (t) - rj(t)|] 



(8) 



and the compression rate £jj = d£ij(i)/dt, where Ri is 
the radius of particle i and fi(t) is its position at time t. 
The expression in square brackets in Eq. ^ may become 
positive during the expansion phase, that is, the (posi- 
tive) dissipative force may overcompensate the (negative) 



III. SIMULTANEOUS CONTACTS 
A. Equations of Motion 

The ICM fails if we take into account the finite dura- 
tion of the collisions. In this case, it may happen that 
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the collision of the particles (process b) starts yet before 
the collision of the lower particle with the floor (process 
a) has terminated. In this case we have a three-particle 
interaction of the floor and both balls which cannot be 
resolved using the concept of the coefficient of restitution. 
In this case, the final velocity of the upper particle must 
be determined by integrating Newton's equation of mo- 
tion for the three-particle system which requires the de- 
tailed knowledge of the interaction forces. Consequently, 
we have to solve the set of Newton's equations 

771121+77113 + ^12-^1)1=0 
m 2 Z 2 + 7772.9 _ F12 = ) 



The latter assumption implies the somewhat unphysi- 
cal fact that the lower side of the large sphere (where it 
contacts the floor) is characterized by a larger dissipative 
constant 701 =/= 712 than its upper side (where it contacts 
the smaller sphere). We will justify this assumption in 
App. [C] where we show that the perhaps more plausi- 
ble assumption 701 = 712, implying eqi 7^ £12, leads to 
qualitatively identical results. 

IV. BASKETBALL - TENNIS BALL PROBLEM 
A. Collision Sequence 



where is the model-specific interaction law between 
particles i and j and the floor is considered as particle 
(with mo — >• 00). 

The failure of the simplifying ICM was discussed in the 
context of the closely related problem of Newton's cradle. 
A simple analysis reveals immediately that the textbook- 
like explanation using isolated collisions is insufficient [S] • 
Instead, the details of the interaction force must be taken 
into account. The explanation of Newton's cradle is far 
from being simple and there is an intensive and contro- 
versial discussion about this seemingly simple classroom 
experiment [91-11 6|. 

The necessity of considering the details of the interac- 
tion force becomes obvious immediately when consider- 
ing colliding rods instead of spheres [T7Ml9| . In fact, the 
investigation of longitudinal waves in colliding bodies and 
the corresponding duration of the collision is a classical 
problem of mechanics, investigated by some of the most 
eminent scientists, such as Poisson [2U], Boltzmann [2T] 
and other important scientists [22H24] 



B. Comparison with the ICM 



Let us assume two vertically aligned balls (the basket- 
ball - tennis ball problem) as sketched in Fig. [I] We 
integrate Newton's equation of motion, Eq. ( |T4| , for this 
system numerically and obtain the forces between the 
bottom and the lower sphere and between both spheres, 
see Fig. [2} For time <t< the lower particle is in 
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In Sec. Ill Bl we conclude the coefficient of restitution 
from the interaction force. Using this result, we can com- 
pute the final relative velocity v' r by means ofEq. O, 
employing the assumption of independent collisions. Al- 
ternatively we can obtain v' r by solving the set of equa- 
tions ( 14 ) numerically. The latter approach does not re- 



FIG. 2: Forces F01 and F12 obtained from solving Eqs. ( 14 1 



quire any assumption on the sequence of the collisions. 
We will see that both results may deviate considerably 
according to rather complex dynamics of the system. 



In order to compare both results by means of Eq. ( 13 ) 



During the contact between the lower particle and the floor 
there occur multiple contacts between the spheres (full line). 
For discussion see the text. Parameters: Ri = 10 cm, R2 = 
lcm, Ah = 0.1mm, z{ 0) = 0.6 m, k = 5.0 • 10 7 N/m, e = 0.7 
(corresponding to hard rubber). 

contact with the floor as indicated by the force Fqi 7^ 0. 
During this interval the balls are in contact repeatedly, 
starting at time (first contact) and ending (last con- 

(e) 



we map the constants of the force law to the coefficient tact ) at time *ia as indicated by F 12 ^ 0. An interesting 
of restitution (fe, 7, m eff ) «-» e. 



We assume that the collisions between the lower sphere 
and the floor (ij — 01) and between the spheres (ij = 12) 
take place at the same coefficient of restitution. Since the 
effective mass enters Eq. (131, for given material stiffness, 
k = const., Eq. (13) then provides a relation between e 



and jij, thus, we can determine 7^ by specifying e as a 
control parameter of the problem. 



detail is the discontuity of F i at t = which is a 
consequence of the force law, Eq. Q: At the instant of 
the contact where £01 ~~ ► 0, the elastically restoring term, 
fc£oij vanishes whereas the (repulsive) dissipative term, 
7^01, has a finite value as soon as the particles get into 
contact. 

The existence of multiple collisions shown in Fig. [2] 
shows that the ICM described in Sec. [n]fails for the cho- 
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sen set of parameters which provokes mainly two ques- 
tions: 

1. How many contacts between the spheres occur and 
how does their number depend on the system pa- 
rameters (Ah, Rx, i?2 7 lij or e respectively)? 

2. If multiple collisions take place, when does the col- 
lision sequence terminate? 

Depending on the system parameters we may obtain 



4*2 — *ai or t\2 > *0i ' therefore, the second question 
must be answered by a definition: The collision-sequence 
terminates at time t = tt (see Fig. [2l when the last con- 



x(e) 



(e) 



.(e) 



, befoi 



tact between the spheres ceases, before the large sphere 
collides with the floor for the second time. 

To answer the first question, we refer to Fig. [3] which 
illustrates the sequence of collisions in dependence of the 
coefficient of restitution e for fixed Ah and i?i/i?2- The 
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FIG. 3: Sequence of collisions for varying coefficient of resti- 
tution. The dashed line shows the end of the contact between 
the lower sphere and the floor, . The fat line corresponds 
to the force drawn in Fig. [2] (Ri — 10 cm, Rz = 1 cm, 



Ah = 0.1 mm, 



,(0) 



0.6 m,k = 5.0- 10 7 N/m) 



value of e was adjusted by varying 7, according to Ei 
(131 while keeping k = 5.0 • I0 7 N/m invariant. Figure p 



should be read horizontally (for fixed value of e): each 
black or grey line marks time intervals when the particles 
are in contact. For elastic balls, e = 1, and the chosen 
parameters there occur 3 collisions. For sufficiently large 
Ah, this number may be unity, that is, the independent- 
collisions condition is fulfilled (see below). Keeping Ah, 
k and the radii R\ and R2 constant and decreasing e, 
the number of contacts increases. This is due to the fact 
that the relative velocity of the balls decreases because 
of inelastic collisions and, thus, the intervals of free flight 
become shorter while the duration of the contacts de- 
pends only weakly on the value of the inelasticity. For 
yet smaller e the relative velocity after the k th contact 



may be small enough such that the lower ball catches up 
with the upper because of its upwards acceleration due 
to its contact with the floor. This effect makes some free- 
flight intervals vanish for decreasing e and, thus, reduces 
the number of contacts. Summarizing, for each set of 
parameters {Ah, R\, R2, k} the number of collisions as 
a function of e is a function with a single maximum. 

For the force law Eq. Q the basketball - tennis ball 
problem may be solved analytically by a piecewise pro- 
cedure, see App. |Aj To check against numerical errors, 
the horizontal gray lines (in between the black lines) show 
the same information as the result of an analytical theory 
which agrees perfectly with the numerical data. 

There is an interesting case when the final velocity of 
the lower ball after losing contact with the floor is only 
slightly larger than the velocity of the upper ball after 
the previous collision. Since both balls move only under 
the action of gravity, the balls may collide an ultimate 
time even after the contact between the lower ball and 
the floor has already finished. These events may be seen 
in Fig. [3] as narrow spikes at e « 0.78, e ~ 0.67, etc. 

The number of contacts of the spheres as a function of 
e and the initial distance Ah is shown in more detail in 
Figure [1] (top) . As explained above for each value of e 
there is an interval for Ah which maximizes the number 
of contacts. 

For the interaction force, Eq. Q the ratio i c .oiAc.i2 of 
the contact duration of the collision lower sphere/ground 
t Ct oi and the contact duration of the collision lower/upper 
sphere t c i 2 increases with m\/m2 (or Ri/R 2 , respec- 
tively). Consequently, the number of contacts increases 
with R1/R2, shown in Fig. [4] (bottom). On the other 
hand, increasing R\/R 2 also increases the initial relative 
velocity between the two spheres and with that the in- 
tervals of free flight, what in turn reduces the possible 
number of contacts. Whereas the effect explained first 
is dominating, the interplay of both effects explaines the 
rather complex behaviour shown in the bottom panel of 
Fig. [4] 

From Figs. [3] and [4] we see that for a vast range of 
parameters the true collision scenario as obtained from 
the integration of Newton's equations of motion deviates 
drastically from the independent-collision scenario out- 
lined in Sec. HU 



B. Effective Coefficient of Restitution 

By solving Newton's equation, we can compute the 
final relative velocity v' r = V2(tf) — v\(tf) which corre- 
sponds to Eq. Q obtained from the ICM. To compare 
both results, we compute v' by integrating Eq. ( Jl4| ) using 
the interaction force Eq. frf\ for a certain set of parame- 
ters {Ah, mi, 77i2, k} and a specified £ = £ spe c (which in 
turn determines 7 via Eq. ( |13[ )). Then, by inverting Eq. 
Q we determine the coefficient of restitution e = e e ff 
which would yield the same final relative velocity for the 
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FIG. 4: Number of contacts between the spheres as a function 
of e and Ah (top) and e and R1/R2 (bottom). The dashed 
lines show the value of Ah and R1/R2 used in Fig. [3] the 
+ symbol shows the parameters used in Fig. [2] The solid 
line in the lower panel indicates R1/R2 = 1. (Parameters: 
k = 5.0 • 10 7 N/m, R 1 = 10cm and R 2 = 1cm (top) and 
k = 5.0 • 10 7 N/m, R 2 = 1cm and Ah = 0.1mm (bottom)) 
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FIG. 5: £ e ff/£s P ec as a function of A/i and the corresponding 
number of contacts (right axis, top). The bottom figure shows 
a magnification of the small Ah range. Parameters: Ri = 
10cm, R 2 = 1cm, k = 5.0 ■ 10 7 N/m, e spcc = 0.9. 



ICM. If £ e ff/£ S pcc ~ 1, both models yield the same re- 
sult, that is, the ICM is an acceptable approximation. 
Otherwise, the ICM fails. 

Consider the dependence of £ e ff/£ S pec on the initial dis- 
tance Ah. For large Ah the lower sphere leaves the floor 
before it contacts the upper one, that is, the ICM holds 
true. Figure [5] (top) shows that £ ff/£ S p C c 1 with in- 
creasing Ah. JVloreover, as expected for £ e ff/£spcc — > 1 
there is only one contact which is a necessary (but not 
sufficient) precondition for independent collisions. Figure 
[5] (bottom) is a magnification of the range of small Ah. 
As discussed before, the number of contacts as a function 
of Ah has a maximum. The oscillations in the number 
of contacts as a function of Ah for very small Aft. cor- 
respond to the spikes shown in Fig. [3] where the lower 
sphere catches up with the upper after the lower sphere 
has already left the ground. 

Similarly, Fig. [6] shows £ c ff/£ S pcc as a function of 
i?i/i?2- Since increasing Ri/R 2 increases the number of 
contacts, £ fr/£ S pcc decreases and thus, as expected, the 




5 10 15 
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FIG. 6: £off/£s P ec as a function of R1/R2 (with R2 = lcm) and 
the corresponding number of contacts (right axis). Further 
parameters: k = 5.0 ■ 10 7 N/m, £ sp0 c = 0.8, Ah = 0.1mm. 



ICM invalidates with increasing R1/R2. 

While for most values of the parameter space 
£ e ff/£ sp ec < I7 there is a small interval of R1/R2 where 
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£cff/£spcc > 1 (inset in Fig, [6]) which is a deviation 
from the ICM too. Here, because of the similar dura- 
tions of the contact upper sphere/lower sphere and lower 
sphere/floor, the lower sphere still being in contact with 
floor pushes the upper one upward. As far as we see, this 
is the only (tiny) effect which allows for e e g > £ spcc . 



V. CONCLUSION 

We considered the motion of two vertically aligned 
spheres which are released to collide with the floor under 
the action of gravity. For the analysis of the dynamics 
of this textbook problem (basketball - tennis ball prob- 
lem) we used two complementary methods. First we de- 
scribed the system exploiting the independent collision 
model (ICM) which assumes instantaneous collisions be- 
tween the spheres and between the lower sphere and the 
floor. The collisions are described by a single number, the 
coefficient of restitution, e, and the duration of the colli- 
sions is neglected. Second, we described the dynamics by 
analytically and numerically solving Newton's equation 
of motion. Here the collisions are characterized by an 
interaction force law, F(£,£). For the case of the linear 
dashpot model used here, the force is a function of the 
elastic and dissipative parameters k and 7. Since there 
is a direct relation between e and {A:, 7}, we can com- 
pare the results of the ICM and the solution of Newton's 
equations. 

We specify two characteristics of the process, a) the 
final relative velocity, v' , between the spheres and b) the 
number of collisions between the spheres during the pro- 
cess. If the approaches were equivalent, we should obtain 
equivalent results for a) and b). 

Obviously, in case of the ICM there is only one colli- 
sion between the spheres and the final result for v' is the 
solution of a textbook problem, Eq. [4] In the article we 
show that Newton's equations yield a different scenario, 
including multiple collisions and the ICM is only valid in 
a certain limit, that is, the ICM fails for a wide range of 
parameters. 

To quantify the deviations, we solve Newton's equa- 
tions with parameters k and 7 that correspond to a cer- 
tain specified coefficient of restitution e = £ spC c- Then we 
compare this value with the effective coefficient of restitu- 
tion, £ c ff, obtained from the final relative velocity as ob- 
tained from Newton's equation. The value e s/e spec = I 
would indicate that both models agree. Our results re- 
veal, however, a dramatic deviation from this ideal be- 
havior. In Figs. [5] and [6] we see that in contradiction to 
the ICM, the ratio £ e s/s spec may adopt any value from 
almost zero up to slightly larger than one, that is, the 
ICM fails dramatically. 

While our subject, the basketball - tennis Ball prob- 
lem, is only a cute but relatively unimportant toy prob- 
lem, our results may have serious consequences for nu- 
merical simulation techniques of granular many-particle 



problems. There exist two established simulation tech- 
niques for the simulation of granular systems, Molecu- 
lar Dynamics (MD) and Event-driven Molecular Dynam- 
ics (EMD). While MD solves Newton's equations of mo- 
tion for all N particles constituting the granular system, 
thus, solves a system of 37V (without rotation) coupled, 
strongly non-linear differential equations, EMD describes 
the dynamics of the ./V-particle system as a sequence of 
pairwise collisions. The latter approach allows for a great 
speedup of the numerical simulation since instead of solv- 
ing computer-time intensive solutions of differential equa- 
tions, we only have to compute postcollisional velocities 
from the precollisional ones, as a function of the coef- 
ficient of restitution for each pair of colliding particles, 
{t7j, v.,-, e} — ¥ {v!,Vj} via a simple propagation function. 
In between the collisions the particles follow simple bal- 
listic trajectories. 

It is obvious, that EMD allows for very efficient sim- 
ulations as compared with MD, in particular for large 
N ~ I0 6 . . . 10 8 , however, this speedup comes for the 
price of the assumption of independent collisions, that 
is, EMD assumes instantaneous collisions neglecting the 
duration of the collisions. While this assumption may be 
justified in a granular gas where the mean free flight time 
is large as compared to the typical duration of collisions, 
it fails for dense systems. Our simple one-dimensional, 
3-particle system shows that the failure may be dramatic. 

For the analytical calculations presented in this arti- 
cle we made two major assumptions whose justification 
might not be obvious beforehand: First we assumed a 
linear-dashpot force, Eq. [7J for the interaction of vis- 
coelastic spheres. This force allows for a simple mapping 
of the constants k and e to the coefficient of restitution 
which is, moreover, a constant in this case. Of course, the 
interaction of spheres is described by a (modified) Hertz 
law which leads to a impact velocity dependent e. We 
could perform the entire calculation presented here also 
for the Hertz law, however, at a much larger mathemat- 
ical effort (see [25] for a similar calculation). We prefer 
here the simplified force and demonstrate in Appendix |B| 
that the Hertz law leads to qualitatively identical results. 

The second simplification concerns the assumption of 
a universal coefficient of restitution for the description of 
the collisions between the particles and between the lower 
particle and the floor. Since the effective mass enters the 
mapping between the force constants and e, the assump- 
tion of a universal e implies that the lower sphere is char- 
acterized by a certain set of parameters {k, 7} when col- 
liding with the floor, but by a different set of parameters 
when colliding with the upper sphere. The alternative 
assumption of invariant material parameters is, perhaps, 
more plausible but leads then to different values of the 
coefficient of restitution for particle-particle and particle- 
floor collisions. While these alternative assumptions lead, 
of course, to different results, in Appendix [C] we demon- 
strate that the qualitative properties of the dynamics are 
the same for both assumptions. 



7 



Appendix A: Analytical description 

For an approximative analytical description, we as- 
sume that the motion of the large ball is not affected by 
the small ball. This adiabatic approximation becomes 
exact for Ri » R 2 and Eqs. (14 1 decouple and may 
be solved piecewise. We obtain four different types of 
collective motion: 

Type A): The balls are isolated from one another and 
from the floor. Here the particles move along ballistic 
trajectories 



zi,free(*) = -gt z /2 + v[ k) t + z\ 
z 2Mc (t) = -gt 2 /2 + v { 2 k} t + z { 2 k) 



(Al) 



In our notation and stand for the positions and 
velocities of the particles at the instant when the system 
enters a type of motion for the fc th time, that is, they are 
initial conditions of the piecewise analytical solution. 

Type B): The lower ball is in contact with the floor while 
the upper one moves freely. Here the upper sphere moves 
along a ballistic trajectory, z 2j { rcc (t) while the lower one 
moves according to a damped harmonic oscillator, 

m\i\ — —m\g + h\{R\ — z%) — "i\Z\ (A2) 

with the solution 



zi,ground(*) = Acos[ujt + p]e + z Unh . (A3) 



where 



p — arctan 



\(z[ k) -z 1Mh )+v[ k) 
(2l,inh ~ z l V 



, _ 4 fc) - gijnh . _ k 1 R 1 -m 1 g 

^ / \ ; ^l,inh — i 

cos(p) fcl 

x - L ; w = A 2 ; uj Q = V fc i/ m i 
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Tj/pe Cj: TTie 6a/fe are in contact and the lower ball con- 
tacts the floor. Here the lower sphere moves due to Eq. 



(A3), disregarding the force resulting from the contact 
with the upper sphere (adiabatic approximation). The 
latter moves like a damped harmonic oscillator in the 
presence of gravity, additionally driven by the motion 
2M,ground(0 of the lower sphere: 



m 2 z 2 



-m 2 g 



k 2 [Ri + R 2 — (z 2 — zi !groun( j)] 

- 72(^2 - h, ground)- (A4) 



The solution z 2 , g round(£) of Eq. (A4) is straightforward 
and similar to Eq. ( |A3[ ) but since it is rather lengthy it 
is not given here. 



Type D): The balls are in contact with each other but not 
with the floor. Here the lower sphere follows a ballistic 
trajectory disregarding the force resulting form the con- 
tact with the upper one (adiabatic approximation) and 
the upper sphere moves as described in type C) but now 
driven by the motion Z\f ree (t) of the lower sphere: 

m 2 z 2 = -m 2 g + k 2 [Ri + R 2 - (z 2 - z li[rcc )} 

- 7 2 (i 2 - il.free) • (A5) 

Again we do not provide the lengthy but straightforward 
solution z 2 ■ 



We keep in mind that ba 



(t) of Eq. (jA5j) here. 

Is i and j (with i 



repre- 



senting the floor) are in contact if the mutual compression 
£ij is positive and the interaction force Fij is repulsive. 
Then we obtain the analytical solution of the problem, 
zi(t) and z 2 (f), from combining the analytical solutions 
of the cases A-D by means of the following scheme: 

1. Type A motion until the lower sphere touches the 
floor at T bcl where (£oi(T b ci) > 0)A( J F 1 O i(T bc i) > 0). 

2. Type B motion until the spheres contact each other 
at T bcu where (6 2 (T bcu ) > 0) A (F 12 (T bcu ) > 0). 

3. Type C motion until the contact between the 
spheres breaks at T GCU where Fi 2 (T ocu ) < 0. 

4. Repeat steps 2 and 3 until the lower sphere leaves 
the floor at T ec i where F i(T ec ;) < 0: 

If Tbcu > T ec i — > Type B motion until T ec i. 
If T ecu > T ec i — > Type C motion until T ec i and then 
Type D motion until the spheres separate at T ccu 
where F 12 (T CCU ) < 0. 

5. Type A motion until 

(a) The lower sphere contacts the floor for the 
second time at T bc i where (£oi(7b c i) > 0) A 
(i 71 oi(T bc i)>0),or 

(b) The spheres touch each other again at T bcu 
where (£i 2 (T bcu ) > 0) A (F 12 (T bcu ) > 0). 

In the first case the collision sequence has termi- 
nated. In the second case: Type D motion until the 
spheres separate at T ccu where Fi 2 (T ccu ) < or 
until the lower sphere contacts the ground at T bc i 
where (£01 (Tbci) > 0) A (F 01 (T bcl ) > 0). 

The described procedure seems to be circumstantial 
but it provides an exact analytical solution of the prob- 
lem in adiabatic approximation. 



Appendix B: Validity of the Simplified Force Model 



The purpose of this Appendix is to demonstrate that 



the analytical and numerical results presented in Sec. IV 
are more than just artifacts of the simplified interaction 
force Eq. 0. The reason of the deviation of the effec- 
tive coefficient of restitution from the specified coefficient 



shown in Figs. [5] and [6] are the described multiple col- 
lisions arising from the finite duration of the collisions. 
Therefore, here we show that multiple collisions also ap- 
pear for the much more realistic interaction force Eq. 
( |B1[ ). To this end we reproduce Fig. [3] where the time 
intervals of particle contacts are shown in dependence of 
the (specified) coefficient of restitution. 



1. Viscoelastic Spheres 



is convenient as it leads to 
Eq. (131, in an elementary 



The force law, Eq. Q 
a coefficient of restitution 
way. However, this force law is a strong simplification. 
Perhaps the simplest particle interaction model which is 
not in conflict with basic mechanics of materials, is the 
contact of viscoelastic spheres, i and j [BJ, given by 



,eff 



with 



3(l-^ 2 ) : 



k = 



P 

„eff' 
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3 pA 



.off' 



(Bl) 



(B2) 



the Young modulus Y, the Poisson ratio v, the effec- 
tive radius = RiRj/ (Ri + Rj) and the effective mass 
rrijj = mivrijl (rrii + rrij). Again we use the mutual com- 
pression and the compression rate to describe the contact 
dynamics (see Eq. The dissipative constant A is a 

function of the material viscosity; see [BJ for details. For 
elastic spheres, A = 0, we recover the classical Hertz 
contact force [25] , 

A necessary prerequisite for deriving Hertz's law of 
contact and its generalization to viscoelastic spheres, Eq. 
(Bl I, is the assumption of small particle deformation, 



that is, the interaction force causes only local displace- 
ments in the region of the contact area. Moreover, the 
impact rate must be small as compared to the speed of 
sound to allow for a quasistatic approximation, see [6j. 
More complex deformations including surface waves and 
oscillations, e.g. [27], are not considered here. Such oscil- 
lations may also give rise to multiple collisions between 
particles and yet more complicated particle-particle in- 
teraction. 

The relation between the coefficient of restitution and 
the parameters of the force law, corresponding to Eq. 
(13), may also be obtained for the case of viscoelastic 

(a sim- 
]), here 



spheres. The calculation is cumbersome [25J [2 
plified version is based on a dimension analysis 
we present only the result: 



£ (v 



OC 

E 

fc=0 



fe/10 



with the initial conditions £jj(0) = 0, ^y(O) = 
with (3 = 7fc -3 / 5 and the pure numbers ho = 1, 



(B3) 

= v and 
hi = 0, 



h 2 = -1.153, h 3 = 0,h A = 0.798, h 5 = 0.267, . . . (see [25J 
for the numerical values hk). Note that in contrast to the 
previous case, Eq. ( 13 ), here the coefficient of restitution 



depends on the impact velocity v. 



2. Basketball - Tennis Ball Problem for 
Viscoelastic Balls 

Just as in Sec. II VI we use the coefficient of restitution e 
to characterize the system's dissipative properties, due to 
the dissipative constant A in the force law, Eq. (Bl ). We 
proceed on the lines of Sec. |IV A[ We specify e, the Young 
modulus Y, the poisson ratio v, the effective radius Rf^ 
and mass m|? and solve Eq. 



(B3) numerically for the 



dissipative parameter A. Additionally, we specify the 
initial velocity Vi„ = 3 m/s since for the viscoelastic force 
law the coefficient of restitution depends on the impact 
As in Sec. 



rate. 



IV 



the assumption of a universal value 
of the coefficient of restitution to describe both particle- 
particle and particle-floor contact results in the rather 
artificial fact that the spheres cannot consist of the same 
material. This assumption is necessary to use the quite 
descriptive quantity e as a control parameter. In App. 
[C]we will release this assumption and show that it does 
not qualitative change the system's behavior. 

Figure [7] shows the sequence of collisions, correspond- 
ing to Fig. [3] for the linear-dashpot force. The figures 
reveal the same structure of the collision scenario, that 



is, the viscoelastic force, Eq. (Bl), leads to qualitatively 



the same results as the linear-dashpot model and, thus, 
justifies application of the simplified force Eq. |7]) in Sec. 




FIG. 7: Same as Fig. [3] but for viscoelastic spheres. Pa- 
rameters: (Ri — 10cm, R2 — 0.25cm, Ah — 0.1mm, 
Y = 5.0 ■ 10 7 N/m 2 , v = 0.45). The black line shows the 
end of the contact between the lower sphere and the floor for 
each e. 
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FIG. 8: Same as Fig. [7] but for particles made from the same 
material, characterized by the dissipative material parameter 
A. (Parameters: Ri = 10cm, R2 = 0.25cm, Ah — 0.1mm, 
Y = 5.0 ■ 10 7 N/m 2 , v = 0.45). 



spheres occur via the same coefficient of restitution e 
which allows to consider e as a control parameter. This 
assumption, however, implies also different dissipative 
constants for the contacts. 



In this Appendix we reproduce Fig. [7J once again with 
the complementary assumption of identical material pa- 
rameters which implies in its turn different coefficients 
of restitution £01 for the particle-floor and £12 for the 
particle-particle contact, see Fig. [5J Thus, we can no 
longer use e as characteristic value. Instead, to char- 
acterize the interaction we use the dissipative material 
parameter A which enters the force law, Eq. (Bl I. 



Appendix C: Validity of the assumption of a 
universal coefficient of restitution 

For the calculations we assumed that the collisions be- 
tween the lower sphere and the floor and between the 



The sequence of collisions shown in Fig. [3] has the same 
structure as for the assumption of a universal coefficient 
of restitution with only minor quantitative differences. 
Still there occur multiple collisions and consequently the 
effects described in Sec. |IV| persist. Hence the assump- 
tion of a universal coefficient of restitution is justified. 
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